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Abstract 

The dynamic charge susceptibility and the optical conductivity are calculated 

in the planar t—J model within the memory function method, working directly 

in terms of Hubbard operators. The density fluctuation spectrum consists of 

a damped sound-like mode for small wave vectors and a broad high energy 

peak (~ t) for large momenta. The study of the optical conductivity shows 

that electron scattering from spin fluctuations leads to the Drude-frequency 

dependent relaxation rate which exhibits a crossover from w 3 / 2 behavior at low 

frequencies (u < 2|/x|), to a linear w-dependence for frequencies larger than 

2\/j,\. Due to the spin-polaron nature of charge carriers, extra absorptions arise 

starting at a frequency ui > J. The obtained results are in a good agreement 

with exact diagonalization studies. 

PACS numbers:71.45.Gm, 73.20.Mf, 78.20.-e 
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I. INTRODUCTION 



Among other unconventional normal state properties of high-T c superconductors, an 
anomalous charge dynamics has also been detected in the optical measurements of the un- 
derdoped samples. Namely, a non-Drude fall-off of the low-frequency absorption indicating 
a linear (^-dependence of the relaxation rate and an anomalous mid-infrared (MIR) band 
with a typical energy ~ 0.1 eV have been observed!!!! 

It is widely believed that unusual properties of the superconducting cuprates are due to 
the strong electron correlations.0 The minimal model to describe correlation effects in the 
cuprates is the t — J model. While a number of analytic works have been carried out to 
investigate spin dynamics within the t — J model, only few of authors have studied charge 
dynamics.H In Refs. HJ5| charge fluctuations have been studied by the slave boson and 
Hubbard operator (HO) formalism within the leading order of 1/N expansion, respectively. 
It was found that the density fluctuations at large momenta show a sharp high-energy 
peak corresponding to the collective mode which reduces to the sound mode in the long- 
wavelength limit .Si Later, the authors of Ref. ^| showed that next-order corrections in the 
1/N expansion lead to the broadening of the high-energy peak due to incoherent motion of 
bear holes. Similar features of the density response have been previously observed in exact 
diagonalization studies of small clusters.! Recently, the charge correlations in the AFM phase 
of the t — J model has been investigated by commulant version of projection technique!. 

In the present paper, we investigate the charge fluctuation spectrum of the t — J model 
in the paramagnetic state with short-range antiferromagnetic (AFM) correlations. We de- 
velop a self-consistent theory for the dynamic charge susceptibility (DCS) by applying the 
memory function method in terms of HO's. The employment of the HO technique has a 
twofold advantage; By using the equations of motion for the HO's we automatically take into 
account scattering of electrons on spin and charge fluctuations originated from the strong 
correlations, as it has first been pointed out by Hubbard.i Moreover, the HO formalism 
allows us to preserve rigorously the local constraint of no double occupancy. 



2 



We calculate the memory function within the mode coupling approximation (MCA) in 
terms of the dressed particle-hole (p-h) and spin fluctuations. Similarly to the nearly anti- 
ferromagnetic Fermi liquid approach,0 we treat fermionic and localized spin excitations as 
independent degrees of freedom. We show that the memory function involves two contribu- 
tions. The first one stems from the hopping term and describes a particle-hole contribution 
from the itinerant hole subsystem. The second one involves scattering processes of electrons 
on charge and spin fluctuations and comes both from kinematics and exchange interactions. 

Further, we perform an analytic analysis of different limiting behavior of DCS to show 
that the essential features observed in the exact diagonalization studies^ can be reproduced 
within the present formalism. We find out that for small q the DCS is mainly governed 
by the sound mode. Although the unrenormalized sound velocity is larger than the Fermi 
velocity, unlike the Fermi liquid theory, the "self-energy" corrections lead to softening of the 
sound. The renormalized sound falls down into the p-h continuum getting a finite damping 
due to the decay into pair excitations. In the short- wavelength limit, density fluctuation 
spectrum mainly consists of a broad high-energy peak. At large enough wave vectors the 
peak is dispersed out of the coherent p-h continuum and broadens due to high energy (~ t) 
transitions involving the incoherent band of one-particle excitations. 

We also discuss the optical conductivity a(uj). For low frequencies we analyze ct{uj) in 
terms of the generalized Drude law. We show that there is a mass enhancement of order 
m*/m ~ 6, due to the electron scattering on spin fluctuations. These scattering processes 
also lead to a frequency-dependent relaxation rate which exhibits a crossover from to 3 / 2 
behavior at low frequencies, 10 < 2|//|, to a linear cu-dependence for lo > 2|/i|. A possible 
origin of the MIR band is also discussed. 

The paper is organized as follows. In the next section, we give the basic definitions and 
sketch the memory function formalism. In Sec. Ill, we employ MCA to calculate a memory 
function. The dynamic charge susceptibility and optical conductivity are discussed in Sees. 
IV and V, respectively. The last section summerizes our main results. 
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II. MODEL AND MEMORY FUNCTION FORMALISM 



The t — J model expressed in terms of HO's, X" 13 = \i, a)(i, (3\, reads as 



H = H t + Hj = -J2 U 3 XfX°° 



+ \Y,JiAxfx™-xrxf}, (i) 

where the indices and a = ±1 correspond to a hole and an electron with spin a/2, 
respectively, tij = t and Jij = J for the nearest-neighbor (n.n.) sites on a planar lattice. 
The HO's can be either Bose-like or Fermi-like and obey the following on-site multiplication 
rules X^X] S = 6 / 3 1 X° iS and the commutation relations 



[S M Xf ± 8 5a Xf) , (2) 



where the upper sign stands for the case when both HO's are Fermi-like, otherwise the lower 
sign should be adopted. In the t — J model only singly occupied sites are retained and the 
completeness relation for the HO's reads as 

x oo + £ Xr = 1 . (3) 

The spin and density operators are expressed by HO's as 

S° = Xf, Sf = WaXr, n l = J2xr. (4) 

The dynamic charge susceptibility -/Vq(u/) is given by a Fourier transformed two-time 
retarded Green function (GF)0 

N^u) = -«n q |n_ q )) w = i / dte^([n q (t),n_ q ]). (5) 

Jo 

To calculate Nq(u>), we employ the memory formalism as discussed in Refs. |l2"Hl4l. First 
we introduce density- density relaxation function 

= ((n q |n_ q )) w = -i / G?te iart (n q (t)|n_ q ), (6) 

J 



where the Kubo-Mori scalar product is defined as 

(A(t),B) = [ dX(A(t-iX)B), (7) 
Jo 

with (3 = 1/T. The DCS iV q (a>) is coupled to the relaxation function <& g (u;) by the equation 

N n (uj) = N (l -LU$ q (u), (8) 

where iV q = iV q (0) is the static susceptibility. 

We introduce the memory function M q (u;) for the relaxation function 3> q (ci>) as 

N 

By adopting the equation of motion method for the relaxation function 3> q (&>) one finds that 
the memory function M q (u;) is given by the irreducible part of the relaxation function for 
"currents"@'@ 

M» = ((j q |j- q ))r- (10) 
The "current" operator j q in the site representation reads as 

ji =fn= -i[n h H\ = -iJ2tij(XfXt ~ H.c). (11) 

j> 

The Heisenberg part of the Hamiltonian ([]]) conserves the local particle number and thus 
gives no contribution to the "current" operator (|TT|) . To treat properly a contribution from 
the Hj term, we go one step further and similarly to Eq.(^) we introduce the memory 
function M q (u;) for the relaxation function for "currents" M2(u) fllQ|) , by the following 
equation 



M Q V) = '" q w , (12) 

qV ; - M q (u;)/m q ' V ; 

where 

m q = -<0"qb'-q))«=0 = iflj'q^-q]) (13) 
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is the first moment of DCS and the memory function M q (u;) is given by the irreducible part 
of the relaxation function for "forces " : 

M q (u,) = ((F q |F_ q ))r (14) 

with 

F q =j q = -i[7 q) tf]. (15) 

Further, to close the system of equations [see Eqs.(|8|),([J), and ([12])], we employ a mode 
coupling approximation for the memory function M q (c<j). 



III. MODE COUPLING APPROXIMATION 

First we express the memory function in terms of the irreducible part of time-dependent 
corelation function for "forces" by means of the fluctuation-dissipation theoremll! 

' dJ eP^-l 



M q (cj) 



2tt u'{u — u' + irj) 

— oo 
oo 

x J die' 1 " '*<F_ q (t)|F q >* 

— oo 



The "force" operator is given by 



(16) 



(17) 



m j,<7<7' 



n: 



in. i, j — tim [ tmj{^i -Aj van 



A,- Orrrjl — Aj Aj -£> m J 



I T ^""O "V^Ocr' rjcr'cr] 



(18) 



where the Bose-like operator 



B° 



A; Oo-o-' — Aj Og-cr' 



-n,- — crS? 



(19) 



describes electron scattering on spin and charge fluctuations. 
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The sum in Eq.(|T7D contains the products of HO's from the same site. As follows, 
such products give no contribution to the memory function, while being decoupled they do 
contribute. That is a result of the complexity of HO's algebra (Q). To show this, let us 
consider the term given by Eq. flT8|) ; since (i,m) and (m,j) are n.n. pairs i^m and m ^ j, 
however i can be equal to j. For the latter case, i = j, the first term in Eq.(|T8D is linear in 
the density operator {X°°X® U = X° a ) and thus gives no contribution to the irreducible part 
of the correlation function for "forces" .0 As for the second term (|18|) , one can easily verify 
that in the case i = j it is canceled out by its counter part from the sum (|l7l). Finally, the 
last term in Eq . fll8|) vanishes for i — j since X[°Bf a = due to the constraint. Therefore 
we have to subtract these terms from the "force" operator. As a result, we come to the 
following expression in the momentum space 

1 



7T7 2^ Ll k,q A k+q A k 

V JV k 



- J2 M K*,p x £+<i-p X k B p °, (2°) 

k,p,cr(j' 

the vertices G k)q and Mk iq , p are given by 

Ck,q = Pk.q - #k,q (21) 
^k,q,p = X][ m k,q,p ~~ m k,q,p 



Mk,q,p=ERq.p-^ ) q,p] ( 2 2) 



where 



( j.\2 2 (i) z 2 tJ r 

2k, q = {zt) 7 k q , m kqp = — -7 P [7k-p,q - 7k,, 



m 



(2) _ ^2 



Zt) [7k,q7k+q-p-7k-p,q7k]- (23) 



In Eqs.(pT|) and (^) g(m) denote g(m) averaged over the Brillouin zone and are given by 

ztJ, 



9k. 



q 



2zt 2 (l - 7q ), m^ ip = — [7k, q - 7k- P , 



= 2^ 2 7p-q,q, (24) 



where 7 k , q = 7k+ q — 7k, 7q = l/2[cos(g x ) + cos(q y )] and z = 4 for a 2-dimensional square 
lattice. This form of the renormalized vertices (|2lD-(p2|) insures that all the operators in the 



products of Eq. ( p0|) are from different sites. Therefore HOs can be simply permuted within 
the decoupling procedure. 

To calculate the irreducible part of the time-dependent corelation function in the right- 
hand side of Eq. (|T6D , we apply the mode-coupling approximation^ in terms of an indepen- 
dent propagation of dressed p-h pairs and charge-spin fluctuations. The proposed approxi- 
mation is defined by the following decoupling of the time-dependent correlation functions 

(X?-*(t)X*'{t)\Xg^X&') ~ 

5a, ( r'4-q,k'(^k q( i )^')( X k <T ( i(; ) X k'+q)> ( 25 ) 

V*k-q-p(*) A k Wtip l*Jl A k'+q-p' A k' -Dp' / 
— ^(j,s'^<j',s^k-q-p,k'5p,-p' 

(^k° q -pW^>(^(0^ +q -P')(< ff (*)<' S )- (26) 

By using the decoupling scheme and the spectral representation for the two-time retarded 
GF's0 we obtain for the memory function 

M q H = -[n(q,w)-n(q,0)], 

UJ 

n(q, W ) =n 1 (q,^) + n 2 (q,^), (27) 



where Ili(q, tu) and IT^q, u) stem from the first and the second term of Eq. (|20|) , respectively. 
Their imaginary parts are given by 

no 

-2tt 



— Z7T f 

n?(q,w) = J du}lU < 



x AklwO^K + ij), (28) 

2 °° 
^(q, w) = M k,q, P / / ^l^-Nowi.wa 

k .P -oo 

X v4 k+q _ p (^ - Wi + W 2 )^k(^2)Xc S (P> W l)> ( 29 ) 

where n Wli£J = n(o»i) - N W:UJUUJ2 = [1 + N(ui) + N(uj — u} 1 )]n U2tW - 0Jl with n{uj) and 

iV(a>) being Fermi and Bose distribution functions, respectively, and 



A h {u>) = --Im((Xj CT | Xf» WJ (30) 

7T 

is a single-particle spectral function which does not depend on spin a in the paramagnetic 
state. We have also introduced the unified spin-charge fluctuation spectrum 

x^(q^) = ^M + x», (3i) 

with 

XM = — «S q |S_ q » tt (32) 

7T 

being the dynamic spin susceptibility. In obtaining Eq. (p9|) we have also used the identity 
((SZ\SZc))u) — 2((S r q|5'iq))aj which holds in the paramagnetic state. 

As it is clear from Eq.(|27|) the memory function involves two contributions; the first 
one ( p8|) stems from the H t term ([I]) and describes the p-h contribution, while the second 
one comes both from H t and Hj parts and involves electron scattering on spin and charge 
fluctuations. Further we assume that the charge carriers are mainly relaxed by the scattering 
on spin fluctuations and retain only the second term XqO^) i n Eq. (|3T|) . 

For further discussion it is more convenient to integrate out the fermionic degrees of 



freedom (k, o^) in Il^q, u>) (|29|) , which results 



1 °° 



xn" (33) 



where we have introduced an effective spectral function 

-2tt 



r 

nj_ p (w) = ^pH^k,q, P y duixn^ 

^ — oo 



X ^k+q-p(^ + ^l)^k(^l), (34) 

for particle-hole excitations coupled to a particular (p, u) state of spin fluctuations. 

To conclude the section, we calculate the first moment of DCS (|T3"D. By performing the 
commutation between the density and "current" (0) operators we readily get 



m q = 4^iVi(l-7 q ), 



(35) 



where 

^4et:w s > (36) 

q 

is the p-h corelation function. 



IV. DYNAMIC CHARGE SUSCEPTIBILITY 



Equations @,(§), and (|T2| ) result in the following form of DCS 

TYl 

N * {U}) = "^-[n(q^)-n(q,0)]/m q -^' ( 37 ) 

where m q is the first moment of DCS given by Eq.(^), and = m q /iV q is a mean field (MF) 
spectrum for the density fluctuations. The memory function formalism does not provide 
itself the static susceptibility iV q . The latter is calculated within the same approximation 
scheme as for the static spin susceptibility^ which results in the following form of the MF 
spectrum 

Ql = 2zH 2 C cl (l - 7q ) (38) 

where 

C q = ±(l-£) + JV 2 -^Vi(l + *7q) (39) 

and the parameters N m are defined by Eq. fl3"6|) . Here we note that the MF spectrum f2 q 
resembles the dispersion of an undamped collective mode in the charge channel found in the 
leading order of 1/N expansion.0 

In a proper analysis, Eqs. (pT|)-(f29"D and (|37| ) should be treated self-consistently with 
the equations for the single-particle spectral function v4k(u;)lil and the spin susceptibility 
This problem, to our knowledge, can be solved only numerically. Here we show, 
that the main features of charge fluctuation spectrum observed in the exact diagonalization 
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studies! can be, at least qualitatively, reproduced in an analytic way based on the physically 
justified anzats for the single-particle spectral function and the dynamic spin susceptibility 
First, we discuss the one-particle spectral function A^iuj). 

Actually, the spectral characteristics of the t — J model have been investigated by various 
analytic and numerical approaches.^ Those results led to the consensus that the single- 
particle spectrum involves a narrow quasiparticle (QP) band of coherent states and a broad 
continuum of incoherent states. The corresponding spectral function can be represented as 

A k (u) = A c ^ h (uj) + A^ c (u), (40) 

where the QP part is given by 

A c k °\uj) = Z k S(uj-e k ), (41) 

with Zy_ and = — \i being the QP weight and dispersion referred to the chemical 
potential /x, respectively While the incoherent part A™(u) is little affected by the doping, 
the coherent band structure strongly depends on the magnetic background. Namely, in the 
low doping regime (ordered phase) the QP dispersion is determined by the hopping within a 
given AFM sublatticelil and the Fermi surface (FS) consists of small hole-pockets centered 
around (±7r/2, ±7r/2). While for a moderate doping (paramagnetic state) the dispersion 
reflects the dominance of n.n. hopping, there exists large electronic FS which encloses a 
fraction of the Brillouin zone equal to the electron concentration n.0 For the latter doping 
regime, the exact diagonalization results are well fitted by a simple tight-binding dispersion 
with some effective hopping amplitude t e g which scales with J ( for 5 = 0.1 and J = OAt 
t e fr = 0.24iS). Hence, in the paramagnetic phase we can put = — zt e fT7k- 

A nearly structureless incoherent part is predominantly distributed below the QP band 
(in the electronic picture) and we approximate A^ c (uj) as follows 

Aj?V) = -Oi-w'MWto: + cu'), (42) 

where Wi nc ~ 5t is the incoherent bandwidth and oj' is measured from the bottom of the QP 
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band. The spectral weights of the incoherent band T and QP are provided by the sum 
rules 



n 



and are given by 



k — oo 

1 ? 

- Y, J da>M<*) = ( 43 ) 



26 r - ^W mc , (44) 



1 + 5' (1-5) 2 

where Z = (Z^) is the averaged QP weight and it coincides with that obtained within the 
Gutzwiller approximation.^ 

With the above form of the spectral function (|40| ) - (|42|) , the equal-time corelation func- 
tions Ni t 2 (^9]) are estimated to be 

2 n 
N x ~—nZ, N 2 ~— (45) 

7T Z ZZ 

Finally, we assume that the spin susceptibility Xq.{ u ) dH) is peaked at the AFM wave 
vector Q = (tt, 7r)S Then, at T = 0, we approximate IT^q, uj) (|33|) as follows 

n^(q,w) ~ J duj 1 fl'^ Q (uj-uj 1 )x"(uj 1 ), (46) 

where 



o 



p 

is the local spin susceptibility. Since the detailed form of x( w ) is n °t essential for our study, 
we use the following relaxation-type susceptibility 

XH = ; * i (48) 

1 — lUJjijjQ 



which seems to be in agreement with the exact diagonalization data.Q^I In Eq. ( f4"8|) ujq oc J 
is the energy scale of spin fluctuations and x is the static susceptibility provided by the sum 
rule 
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3 

duj X {uj) = -n, (49) 



with the high-energy cutoff u c = 2J. Equations (HHf) and (Ef9|) result in 



3n 

X = 2^ ln[l + (o; c /a;o) 2 ]- (5 ° } 



A. Long wavelength limit 

First, we discuss a small q, uo limit of DCS. In the long wave-length limit, q — > 0, the MF 



spectrum fl38|) reduces to the sound mode f2 q = v s q with the sound velocity v B = ztJCo/2 



larger than the Fermi velocity vp(jj,) = zt c sy (1 — /i 2 )/2 ( ji = \fi\/zt c s ~ 7r<5/4 for small 5). 
At small q the vertex functions (|21"|)-(|24|) are given by 



G k , q = 2t 2 [2^(qV k7k ) 2 - l}q 2 , (51) 

M kiQiP = z(z - l)t J(qV k 7k)g. (52) 

with q = q/g. 

First we consider the real part of "self-energy" Il(q, uo) (]27|)-(|2"9"|). Since for small q, the 



vertex functions Gk,q ~ <? 2 and M ki Q iP ~ q, we approximate n'(q, uo) ~ II 2 (q, to) to keep the 
contributions leading in small q. Moreover, for small u we expand lT 2 (q, uo) as 

n 2 (q,^)~n 2 (q,0)-a qW 2 , (53) 

where a q > and is given by 



^ _ ld 2 IT 2 (q,cj 



1 ^'«|^, (54) 



duo 2 



=0 7T J Co> 



Since n 2 (q, to) d46|) is an odd function of u;, there is no term linear in to in the expansion 



(p3j). Equations (p7j ) and (|53|) result in the following form of DCS for small q,uo 



w - r'"^ (1 + A V ■ (55) 

cj^ — v jq 1 + 2zu;r q 



where 
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v s and T q are the renormalized sound velocity and sound damping, respectively. 

First, we consider renormalization factor A. The estimation of A is interesting by itself, 
since it represents the electron mass enhancement factor [see Sec. V] and can be related to 
experiment.^ Approximation (|40 ) for the one-particle spectral function A^uj) leads to two 
different contributions to the effective spectral function for p-h excitations W'q(uj) ([34]). 
The first one n q _q(u;) c ~ c is due to the transitions within the QP band and the remaining 
part Ilq_Q(u;) 1 " c is provided by the incoherent-coherent transitions. 

First, considering n q _ p (co>) c ~ c ([51]), ( p)ID we come to the following expression for small q 

nUH M ~? 2 E(fe) 2 
iV k 

x Nek) - n(e k + u)]6(u - e k+q _ Q + e k ), (57) 

with A = [z(z — l)tJ] 2 . As it has been previously discussed by several authors.0 the p- 
h spectral function for tight-binding electrons exhibits a crossover at frequency uj = 2\/j,\. 
Namely, for uj < 2\fi\ it is peaked at the incommensurate wave vectors Q* = (tt±5*, tt), (it, 7r± 
5*) where the displacement 5* for small /x is given by 5* ~ f^/t c g, whereas for uj > 2|/x| the 
p-h spectral function gets its maximum value at AFM wave vector Q and follows the nested 
Fermi liquid scalingS Hence we consider the cases uj < 2\fi\ and uj > 2\/jl\ separately. 

In the case uj < 2\fi\ we put e k+q _Q = e k _Q it in Eq . (|57D and for 5* <^ 1 expand e^Q, ~ 
— e k — 5*de] Si /dk x . Moreover, since the dominant contribution to the integral in Eq.(|54]) comes 
from uj ~ 0, at T = we approximate n(e k ) — n(e k+w ) ~ c<j5(e k ) and, as a result, we get 

n q - Q M C " C = ~^ 2 MZ \ 2 uJh{uj)6(2\v\ - uj), (58) 

where = 2zt cS is the coherent bandwidth, and 

2 sin 2 k* 



h(uj) 



7r 2 1 cos k* sin k* \ 



with sink* = 1 — uj/{2\h\) and cos/c* = 2jl — cosk*. For uj 2\\i\ we have h{uj) 
2^Jcj/\fi\/TT 2 and thus 
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— AM 7 2 



^coh V ^ 




Then from Eqs. fl46|) and (|59|) for the velocity renormalization factor A (^g) we get 



V2AZ 



where the QP weight Z and n.n p-h correlator Ni are given by Eqs. (|44j) and (|45|). For the 
actual values of the parameters J = OAt, t e fr = 0.24t, and 5 = 0.1 we obtain A 1 ~ c ~ 4. 

Next we consider the case u > 2\/j,\. In this case the p-h spectral function is peaked at 
the AFM wave vector Q and from Eq. (|57|) we come to 



ICqM C - C = -^q 2 me{u - 2/^(2 - Co), (61) 

Wcoh 



where u = 2uj/W co ^ and 



I(u)) = _L [u + E(u_/io + ) - 2cDK(w_/w+)] (62) 

with C)± — 2 ± u), fT(x) and i?(x) are the complete elliptic integrals of the first and second 
kind, respectively. The function I(x) (|62| ) is normalized to 1/4 in the interval < x < 2 
and well approximated by the following linear dependence I(x) = (2 — x)/8 and we get 

' coh 

/(w) = (iy coh - w)0(u; - 2|^|)6>(W coh - lu). (63) 
That results in the following contribution to A 

H w*E 7 ' (64) 



with 



OO X ~ 

I = [ [ dxdy nX ~ y) x{y)0{2 - y) (65) 



o o 



where the dimensionless functions f(x) and x(y) stand for f(uj) ( |63| ) and spin susceptibility 
x(u;) ( P^ ) measured in units of J. For J = 0.4t and 5 = 0.1 we calculated integral 
numerically and found / = 0.5. Then from Eq.(|64|) we estimate A|~ c — 0.8. 
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As for the remaining contribution H'^_q(uY' c , with the help of Eqs.(p4[),(^2|), and (|52| 
at T = we obtain 

-nAZ o 



AzT 



-q 



x Q(2uj - W coh + 2\ f j,\)e(W tot - u), 



(66) 



where W to t = W co h + Wine is the total bandwidth. From Eq.fl5E|) and by using the sum rule 
(£l9|), the upper value of the incoherent-coherent contributions to fl^q, ujY~ c is estimated as 



nZ( q , w )^-nfi;_ Q (w) 



(67) 



That results in 



A' 



3n(z - l) 2 ZtJ 2 



16riW ( 2 



coh 



4- 



W 2 



coh 



0.1. 



(68) 



A small value of A i_c ~ 0.1 in comparison with A l_c ~ 4 is due to a large threshold 
energy {W co ^/2 — — t) ( |66D for incoherent-coherent transitions which are important only 
in describing high-energy density fluctuations. 

Finally, by summing all three contributions (0), (|66"D , and (|70|) for the velocity renor- 
malization factor we get A ~ 5. Here we notice, that in the Fermi liquid theory, for v s > vp 
"self-energy" corrections stiffen the sonnd.i Contrary to this, , in the present case the sonnd 
softens. That is due to the scattering on spin fluctuations given by the n 2 (q, uj) term (p6|). 
One can easily verify that the renormalized sound velocity (^Bj) gets smaller than the Fermi 
one v s = v s / a/1 + A < vp. It falls down into a particle-hole continuum getting finite damping 
due to the decay into particle-hole pairs.il This process is described by n^q, uj) (^) . The 
latter in the small q, uj limit reads as [see Eqs. (|28|) and (|5T|) 1 

-87rf 4 Z 2 



n;'(q,^ 



N 



-ujq^iM^^-lf 
k 



x 5(e k )8(uj/q - qv k ) 



(69) 



where Vk = Vk^k is the QP velocity. The integration over k in Eq.(|69"D results in 



n;'(q^) 



-V2t 4 Z 2 
nt 2 s \ sinful 



2oj 
av Fi 



ujq 3 , 



(70) 
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where cos 6 = \^i 2 + 2cu 2 /g 2 ]/fp — 1. This gives the following form of the sound damping 

r q ~(3v s q (3= (71) 

for actual values of the parameters /3 < 1, and thus, in accordance with Ref. [], one obtains 
that the sound damping is only numerically smaller than its energy. 



B. Short wavelength limit 

At large momenta the main spectral weight of density fluctuations is located at high en- 
ergies, (~ t), near the MF spectrum f2 q Q58|). For instance, at q = Q we have Qq ~ zt while 
in the exact diagonalization studies^ the peak is observed at to ~ 6t. However, considering 
" self-energy" corrections and noting that W (u) falls off as 1/uj 2 at large frequencies, from 



Eq.d37|) we obtain the renormalized spectrum as f2 q ~ y^ q — n q (0)/m q . One can easily 
show that II q (0) < and hence the spectrum is shifted to higher energies. At large momenta 
the peak is dispersed out of the coherent p-h continuum and its broadening is only due to 
high-energy transitions involving the incoherent band. Since the latter has been neglected 
in Ref. |5], the authors observed an infinitely sharp peak. However, as it follows, the damping 
of the high energy mode is comparable to its energy. Near the pole f2 q we estimate the 
damping as 

r q = — (72) 

2m q H q 

where n"(q, u) = n"(q, u)Y~ c + ^'(q, u>y~ c describes incoherent-coherent transitions. From 
Eqs.(|28|) and (f42"D for q = Q we obtain 

n ; (Q , w) .-^zM^!£ (73) 



The second contribution n^q, u) 1 ' from Eqs.fl34|), (f42|), and ( |6T[) is estimated as 



n;(Q,»)^^. (74) 



Equations ([T§,([73|), and ([Tj) result in 
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f , = 2(25 + 24 »>^ Z ~ 3f. (75) 

Thus, the peak gets rather broad in accordance with the exact diagonalization results.0 

For large momenta but low energies, charge excitation spectrum should show some low 
energy structure related to the contribution from the p-h continuum to N^(u). Since f2 q is 
larger in (£, £) direction than in (£,0), the low energy structure should be less pronounced 
in the latter case. The same anisotropy has been observed in the exact diagonalization 
studies.^ 

V. OPTICAL CONDUCTIVITY 

In this section, we discuss the optical conductivity <j(oj). In the linear response theory 
of Kubol!, the frequency-dependent conductivity is given by the relaxation function for 
currents 

o"«b(w) = — {{J x \Jx))l>- (76) 

By using the continuity equation and equation of motion for the GF's, one can easily relate 
the longitudinal conductivity to the dynamic charge susceptibility (El) 



ie 2 ujN^uj) 

a xx (u) = — — Inn - — , 77 

V 9^0 q l 



where q = q x . From Eqs. (|37D and (|77|) we express conductivity in terms of the memory 
function 

a ^ u) = V^M7) (78) 

where D = limm q /g 2 = ztNi is the Drude weight which is given by one-half the averaged 
kinetic energy D = —(H t )/2, m q and N± are defined by Eqs. ( pop and (p6[), respectively. 
The memory function M(u) reads as 

mm = SM^M, (79) 



where 



(80) 



with II(q, uj) defined by Eq. Q2"7|). Since for small q, Il 1 (q, uj) ~ g 4 and n 2 (q, uj) ~ g 2 , only 
the second one contributes to YI{uj). The latter is given by Eqs.(^)and ( p£3| ) at g = with 
^k,q, P replaced by the transport vertex given by 

M k p = lim Mk ; q,p = [tk-p^k - ^k^k-p - 2tw p ] 

<2^0 g z 

+ ^ + </ P ]k- P -^k], (8i) 

where tk = ztj^, <^k = £^7k> and t>k = dt^/dk x . 

We rewrite conductivity (|78|) in the form of the generalized Drude law as follows 

/ x e 2 D(u) 

V ±/t{uj) — iuj 

where an effective Drude weight and the relaxation time are given by 

~ D 1 1 

M " TTaR f(u) ~ r(w)(l + AH)' (83) 

with 

UJ r { UJ ) 

and 1 + \(uj) is the interaction-induced optical mass enhancement factor. The latter in the 
static limit is calculated in the preceding section and is estimated to be of order 6. That 
is in a good agreement with the optical measurement dataB Optical conductivity of the 
t — J model within the present formalism was studied by one of us (N.P.) in Ref. |14], where 
temperature and frequency dependence of cf(uj) where discussed. Here we mainly focus on 
the analysis of the low-frequency behavior of the relaxation rate 

r(u) = _L = -5!M. (85) 

t(uj) uj 
Following Sec. IV we approximate n"(k;) as 
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, n'A (u -LUl) to < 2\a\ 
U"(u)^ du lX "(u; l )x{ QA m (86) 



IIq(u-ui) lo > 2|//| 

In the case u < 2|/i|, the effective spectral function of p-h excitations IIq^u;) (|59|) is given 
by 

nw * (87) 

We remark the square-root behavior of the p-h spectral function U.Q t {u) ~ y/u instead of 
the conventional linear cj-dependenceB This behavior results in the square-root singularity 
of the structure factor that is known as 2kp anomaly familiar for the electron system in low 
dimensions.il It also leads to the deviation from the conventional square law resulting in the 
following form of the relaxation rate 

r(w) ~ ^^cj 3/2 . (88) 



15n DWl h ^J\fi\uj 

Here we note that the o; 3 / 2 -law of inverse life time for electron states near the saddle points 



has been obtained in Refs. In the present case, the Van Hove singularity plays no role. 

The obtained co> 3 / 2 -dependence of the relaxation time is rather due to the coexistence of the 
peak in the spin fluctuation spectrum and the 2kp "anomaly" in the p-h spectral function 
at q ~ Q. We notice that the former one, 2kp "anomaly" , is not related to the FS topology 
and is inherent in a low-dimensional electron system. 

Now we consider the region u > 2|/x|. In this case the p-h spectral function is peaked 
at the AFM wave vector and is almost (^-independent for low frequencies u <C Wcoh (|63|)- 
which results in 

Unlike the previous case, now the electron band structure is mainly responsible for obtained 
behavior. Of course, the AFM character of spin fluctuations favors the scatter process with 
momentum transfer Q and, thus, enhances its contribution to the relaxation rate. 
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To summarize the low-energy behavior of optical conductivity, we have shown that the 
relaxation rate due to the electron scattering on spin fluctuations exhibits the crossover from 
the T(u) ~ cu 3 / 2 behavior at low frequencies u <C 2|/i| to a linear w-dependence at to > 2|/x|. 

Now we discuss the conductivity at intermediate frequencies. The exact diagonalization 
studies of the latter quantity have suggested a possible explanation of the MIR absorption 
within the one-band modelJilil For instance, as it has been observed in Ref. |32|, the finite 
frequency part of a{oS) is dominated by a single excitation which scales with J in the under- 
doped regime.il The origin of this excitation was ascribed to transitions in which internal 
degrees of freedom of the spin-bag QP are excited. The presence of extra absorption rang- 
ing from MIR frequency to ~ 1 eV was also observed in the 1/N expansion study of the 
t — J model.il The authors of Ref. [33] interpreted this feature as being due to the incoherent 
motion of charge carriers. Since the existence of broad incoherent band in the density of 
states of charge carriers is due to internal degrees of freedom of the spin-bag QP it follows 
that the underling physics of both these points are the same. Bellow we also support this 
explanation of the MIR band. 

Actually, with increasing energy an extra channel of optical transitions opens. These 
are the transitions which involve an incoherent band of the single-particle spectral function. 
As we have already discussed, the incoherent-coherent transitions are characterized by the 
energy scale A = W co h — 2\fi\ being a threshold energy for creating "particle-hole" pairs 
with a "hole" in the incoherent band. Due to this extra channel at uo > A, the real part of 
ct{uj) starts to increase. Since cr{ui) vanishes in the limit to — > oo, there should be a peak 
in conductivity at energies of order A. Since the coherent bandwidth (and hence A) 
scales with J it follows that the typical energy of the peak is also J that coincides with the 



energy scale of fine structure in high-energy absorption found in Ref. |32. 

To conclude this section, we discuss the doping dependence of the Drude and regular 
parts of the optical conductivity. The doping dependence of the renormalized Drude weight 
D fl55| ) is due to the ^-dependence of both kinetic energy and mass enhancement factor. 
For small S the former one scales as | < H t > | ~ ZN\ ~ 5 and increases with doping. 
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While the mass enhancement factor, 1 + A, that is due to the electron scattering on spin 
fluctuations gets smaller upon doping since the system moves away from the AFM phase 
boundary. Combining both points we conclude that the Drude weight increases faster than 
the number of doped holes. Since the sum rule for the optical conductivity has to be 
fulfilled we expect transfer of the spectral weight from the regular to the Drude part of the 
conductivity. The loosing of the spectral weight by the mid infrared band is provided by the 
fact that the normalized rate of the incoherent- coherent transitions U'Kcu) 1 ' / D ~ Z/DT ~ 
(1 — 5) decreases upon doping. This qualitative picture of doping dependence of <j(uj) is 
in general agreement with finite clusters calculations. To give some quantitative estimates 
on evaluation of <t(uj) with doping one needs a more detailed information concerning 5- 
dependence of one-particle spectral characteristics as well as spin fluctuation spectrum used 
as the inputs in our theory. However, the presently available numerical data on the 5- 
dependence of these quantities is far from convincing. 



VI. CONCLUSION 

To summarize, we have developed a self-consistent theory both for the dynamic charge 
susceptibility and the optical conductivity within the memory function formalism in terms 
of the Hubbard operators. In this framework the charge susceptibility has been expressed 
via the memory function M q (u;), that is given by fluctuating force correlation function. 
The essential approximation of the preceding theory is the mode-coupling approximation. 
Within the MCA the memory function has been factorized in terms of single-particle GF 
and spin and charge susceptibilities. As a result a type of the "golden-rule formula" has 
been obtained: the density fluctuation created at some place and time will decay either into 
p-h pair excitations and/or p-h pairs and spin (charge) density fluctuations if the system 
evolves into the future. Though, at present time the precise range of validity of MCA can 
not be given the following general arguments can be used to justify it. First, being a self- 
consistent scheme of calculations MCA does not use any unphysical for strongly correlated 
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systems the zero-order GF and can be viewed as an analog of the non-crossing or the self- 
consistent Born approximation (SCBA) for the single-particle GF. However in the SCBA 
vertex corrections are neglected and only the skeleton loop diagrams are taken into account. 
While MCA correctly reproduces the back-scattering term in the impurity problem^ and 
thus partially accounts for vertex corrections [see Eqs. ([H]) and ( P^ ) for the present case]. 
Next, this approach allows to deal with constrained electron operators avoiding an auxiliary 
field representation that rises the problem of local constraint. Moreover, within our decou- 
pling scheme the operators from the same site is never decoupled and thus the strong local 
correlations are retained. Finally, it is expected that there is no divergence in renormalized 
vertexes signaling the existence of a critical low energy mode and the vertex renormaliza- 
tions that is not accounted by our theory can change only our numerical estimates but not 
the qualitative behavior of the susceptibility. 

Our main findings are as follows. We have shown that in the long-wavelength limit 



the charge fluctuation spectrum is mainly governed by the sound mode flog). Although 
unrenormalized sound velocity is larger than the Fermi velocity, the "self-energy" corrections 
soften the sound. Sound falls down into the the particle-hole continuum and thus acquires 
a finite damping due to the decay into pair excitations. The sound damping flTTJ ) is only 
numerically smaller than its energy and hence there is no well-defined sound mode. 

At large momenta the density fluctuation spectrum mainly consists of a broad high- 
energy peak which nearly follows the MF dispersion (^). At wave vectors large enough 
the peak is dispersed out of the coherent particle-hole continuum and its broadening ([72] ) is 
due to the high energy ~ t transitions involving the incoherent band of the single-particle 
excitations. 

We have also discussed the optical conductivity. At low frequencies we have analyzed 
cr(u) in terms of the generalized Drude law. We have shown that there is a large mass 
enhancement of order m*/m ~ 6, due to the electron scattering on spin fluctuations. This 
scattering process also leads to the non-Drude fall-off of the low energy part of o'iuo). Namely, 
the relaxation rate shows a power law (^-dependence with the exponent 3/2 at low frequencies 

23 



u) < 2\fi\ and it is linear in u at frequencies u > 2|//|. As for the intermediate frequency 
conductivity, we have pointed out the existence of a characteristic energy A (of order J) 
above which an extra channel of the optical transitions opens. These are the transitions in 
which particle-hole pairs with a "hole" in the incoherent band are excited and they might 
be responsible for the experimentally observed MIR absorption. 

The obtained results are in good agreement with the exact diagonalization studies of 
small clusters.0 

In this paper we have focused on the moderate doping regime, when the Fermi surface is 
large and there exist a strong short-range AFM correlations in the system. Our assumption 
about the large electronic FS gives a lower boundary S c of the range of validity of the present 
study with respect to the hole doping (<5 C being the hole concentration when the topology 
of the FS changes, i.e. when a transition from large electronic to the hole pocket-like FS 
takes place when doping becomes smaller than <5 c )il. To discuss the low doping regime 
6 < 5 C one has to take into account the changes in the QP spectra. Here, we only point 
out, that for 5 < 5 C the new low energy scale vp ~ appears in the charge dynamics. 
Due to this and the different topology of FS one can expect the different behavior of the 
low energy charge fluctuations that are driven by FS related particle-hole excitations. The 
existence of the upper boundary is due to the assumption on the presence of strong short 
range AFM correlations. When these correlations are already destroyed the assumption 
that the scattering on spin-fluctuations gives the main input to the charge relaxation is not 
valid anymore and in order to relate the theory with experiment one has to include also the 
scattering on other bosonic degrees of freedom. However, this doping regime is beyond the 
scope of the present consideration. 
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